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Abstract 

We propose a stochastic dynamic model of migration and economic aggregation in a system 
of employed (immobile) and unemployed (mobile) agents which respond to local wage gra- 
dients. Dependent on the local economic situation, described by a production function which 
includes cooperative effects, employed agents can become unemployed and vice versa. The spatio- 
temporal distribution of employed and unemployed agents is investigated both analytically and 
by means of stochastic computer simulations. We find the establishment of distinct economic 
centers out of a random initial distribution. The evolution of these centers occurs in two dif- 
ferent stages: (i) small economic centers are formed based on the positive feedback of mutual 
stimulation/cooperation among the agents, (ii) some of the small centers grow at the expense 
of others, which finally leads to the concentration of the labor force in different extended eco- 
nomic regions. This crossover to large-scale production is accompanied by an increase in the 
unemployment rate. We observe a stable coexistence between these regions, although they exist 
in an internal quasistationary non-equilibrium state and still follow a stochastic eigendynamics. 
KEYWORDS: agents, aggregation, economic geography, stochastic dynamics 



1 Introduction 

In recent years, there has been an increasing interest to link the discussion about complex phe- 
nomena in natural science, such as physics or biology, in particular to those in life science, such 
as sociology, economy or regional planning. One of the challenges is to reveal cross-links between 
the dynamic models used in the specific fields, in order to find out pieces for a common theory of 
self-organization and evolution of complexity. 
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With respect to economic and urban problems, this perspective seems rather new. In his essay 
about self-organization in economics, Krugman (1996b) states: "so far this movement has largely 
passed economics by." However, recent years have also seen numerous attempts for developing a 
new, evolutionary view of economics (Anderson et al., 1988, Wei-Bin, 1991, Arthur et al., 1997, 
Silverberg, 1997, Schweitzer and Silverberg, 1998) 

In a world of fast growing population and intercontinental trade relations, the local emergence of 
new economic centers on one hand, and the global competition of existing economic centers on the 
other hand has a major impact on the socio-economic and political stability of our future world. 
This problem also confronts the natural sciences on their way to a new understanding of complex 
phenomena. This paper aims to contribute to this discussion by adapting a model of interactive 
structure formation, which has already proved its versatility in a variety of applications, to the 
problem of economic agglomeration. 

Spatio-temporal pattern formation in urban and economic systems has been investigated for quite 
a long time (Henderson, 1988, Fujita, 1989, Puu, 1993, Krugman 1991, 1996a). There is also a 
well known German tradition in location theory, associated with names like von Thiinen, Weber, 
Christaller (1933) and Losch (1940). Recent approaches (sparsely recognized especially by American 
economists) tackle the problem of settlement formation, location theory and city size distribution 
within a stochastic dynamical theory, based on master equations and economically motivated utility 
functions (Dendrinos and Haag, 1984, Weidlich and Munz, 1990 a,b, Haag et al, 1992, 1994, 
Weidlich, 1991, 1997). 

Different from these investigations, the current paper focuses on an agent based dynamics of eco- 
nomic concentration. Agent models, which have originally been developed in the Artificial Life 
community (Maes, 1990, Steels, 1995) recently turned out to be a suitable tool for describing eco- 
nomic interaction (Anderson et al, 1988, Marimon et al, 1990, Holland and Miller, 1991, Lane, 
1992, Arthur, 1993, Kirman, 1993, Epstein and Axtell, 1996, Arthur et al, 1997). 

The rational agent model, one of the standard paradigms of neoclassical economic theroy (Silverberg 
and Verspagen, 1994), is based, among others, on the assumption of the agent's complete knowledge 
of all possible actions and their outcomes or a known probability distribution over outcomes, and 
the common knowledge assumption, i.e. that the agent knows that all other agents know exactly 
what he knows and are equally rational. 

In this particular form, the rational agent is just one example of a complex agent, which can be 
regarded as an autonomous entity with either knowledge based or behavior based rules (Maes, 
1990), performing complex actions, such as BDI (belief-desire-interactions) (Miiller et al, 1997). 
The complex agent i.e. is capable of specialization, learning, genetic evolution, etc. However, the 
freedom to define rules and interactions for the agents, could very soon turn out to be a pitfall, 
because of the combinatoric explosion of the state space: for 1000 Agents with 10 rules, the state 
space contains 10 1000 possibilities, hence, almost every desirable result could be produced from such 
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a simulation model. 

The alternative to the complex agent, promoted in this paper, could be the minimalistic agent, 
which acts on the possible simplest set of rules, without deliberative actions. Instead of specializa- 
tion, the minimalistic agent model is based on a large number of "identical" agents, and the focus 
is mainly on cooperative interaction instead of autonomous action. 

A version of the minimalistic agent model, which allows to apply the advanced methods developed 
in statistical physics and stochastic theory, is based on active Brownian particles (Steuernagel et 
al, 1994, Schimansky-Geier et al, 1995, 1997, Schweitzer, 1997, Ebeling et al, 1998). Generally, 
active Brownian particles are Brownian particles with internal degrees of freedom. As a specific 
action, the active Brownian particles (or active walkers, within a discrete approximation) are able 
to generate a self-consistent field which in turn influences their further movement and behavior. 
This non-linear feedback between the particles and the field generated by themselves results in an 
interactive structure formation process on the macroscopic level. Hence, these models have been 
used to simulate a broad variety of pattern formations in complex systems, ranging from physical 
to biological and social systems (Lam and Pochy, 1993, Schweitzer and Schimansky-Geier, 1994, 
Ben- Jacob et al, 1995, Lam, 1995, Schweitzer et al, 1997, Helbing et al, 1997, Schweitzer and 
Steinbrink, 1997). 

In the following, the active Brownian particles are regarded as economic agents with two different 
internal states and a very simple behavior: employed agents, which are immobile, generate a wage 
field (as the result of their work), while unemployed agents migrate guided by local gradients in 
the wage field, i.e. they try to move toward those regions in their vicinity with a high productivity 
and a high value of the wage field. Further unemployed agents can be employed, dependent on the 
local economic situation, but employed agents also can be unemployed. 

In Sect. 2 we present our stochastic approach to the problem of employment and migration. In Sect. 
3, the general model is used to derive a special dynamics recently applied by Krugman (1992, 1996b) 
for economic aggregation. In Sect. 4, we present the economic assumptions involved in our model, 
while Sect. 5 describes the results of computer simulations, which show economic aggregation on 
two time scales. Sect. 6 gives some conclusions and an outlook of possible extension of the model 
presented. 

2 Dynamic Model of Migration and Employment 

Let us consider a two-dimensional system with i = 1, N economic agents, which are represented 
by active Brownian particles. These particles are characterized by two variables: their current 
location, given by the space coordinate j-j, and an internal state, 0j, which could be either one or 
zero: € {0, 1}. Active particles with the internal state 9 = are considered employed agents, Cq, 
active particles with 9 = 1 are considered unemployed agents, C\. With a certain rate, k + (hiring 
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rate), an unemployed agent becomes employed, while with a rate k (firing rate) an employed agent 
becomes unemployed, which can be expressed by the symbolic reaction: 

fe- 

Co ~ Ci (1) 

k+ 

Employed agents are considered immobile, while unemployed agents are able to migrate. The move- 
ment of a migrant may depend both on erratic circumstances and on deterministic forces which 
attract him to a certain place. Within a stochastic approach, this movement can be described by 
the following overdamped Langevin equation: 

d ^ = f( n ) + V2em (2) 

f(ri) describes the local value of a deterministic force, which influences the motion of the agent. 
We note here, that the agent is not subject to long-range forces, which may attract him over large 
distances, but only to a local force. That means the migrant does not count on global information 
which guide his movement, but responds only to local information, which will be specified later. The 
second term in eq. (2) describes random influences on the movement of the individuals modelled 
by a stochastic force £(i), which is assumed to be Gaussian white noise: 

&(*)> = o ; (z i (mj(t')) = s ij 5(t-n (3) 

e is a measure of the strength of the stochastic force. As eq. (2) indicates the unemployed agent will 
move in a very predictable way if the guiding force is large and the stochastic influence is small, 
and he will act rather randomly in the opposite case. 

The current state of the agent community can be described by the canonical iV-particle distribution 
function P[9\,r\, On, vn, i), which gives the probability to find the N agents with the internal 
states 0\ , . . . , On in the vicinity of r\ rjv at time t. The change of this probability can be described 
by a multi-variate master equation (Feistel and Ebeling, 1989), which considers both changes in the 
internal states and migration of the agents. Here, we restrict the discussion to the spatio-temporal 
densities ng(r,t) of unemployed and employed agents, which can be formally derived from the 
iV-particie distribution function: 

N 

n e (r,t) = y^Sgfc S(r - ri) P(0 1 ,r 1 ...,0 N ,r N ,t) dr 1 ...dr N (4) 
i=i jA 

5g t e i is the Kronecker Delta function for discrete variables, which is 1 only for = 0{ and otherwise, 
while <5(r — rj) is Dirac's Delta function used for continuous variables, and A is the system size 
(surface area). With £ {0, 1}, we obtain the spatio-temporal density of the employed agents, 
no(r,t), and of the unemployed agents, n\(r, t). For simplicity, we omit the index 0, by defining: 

n (r,t) = l(r,t) ; ni(r,t) = n(r,t) (5) 
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We only assume that the total number of agents is constant, while the density of employed and 
unemployed agents can change in space and time: 

N (t)= [ l(r,t)dr; N x (t) = f n(r,t)dr; N (t) + N^t) = N (6) 
J A J A 

The density of the employed agents can be changed only by local "hiring" and "firing" processes, 

which with respect to eq. (1), can be described by the reaction equation: 

^l(r,t) = k + n(r,t) - k~ l(r,t) (7) 

The density of the unemployed agents can be changed by two processes, (i) migration and (ii) 
hiring of unemployed and firing of employed agents. With respect to eq. (2), which describes the 
movement, we can derive the following Fokker-Planck equation: 

d d d 2 

— n(r,t)= - — f{r,t)n{r,t) + D n -^n(r,t) 

- k + n(r,t) + k~ l(r,t) (8) 

The first term of the r.h.s. of eq. (8) describes the change of the local density due to the force 
f(r,t), the second term describes the migration of the unemployed agents in terms of a diffusion 
process with D n = e beeing the diffusion coefficient, the third and the fourth term describe local 
changes of the density due to "hiring" or "firing" of agents. 

By now, we have a complete dynamic model which describes local changes of employment and un- 
employment, as well as the migration of unemployed agents. However, so far some of the important 
features of this dynamic model are not specified, namely (i) the deterministic influences on a single 
migrant, expressed by f(r,t), (ii) the hiring and the firing rates, k + , k~, which locally change the 
employment density. These variables depend of course on local economic conditions, hence we will 
need additional economic assumptions. 



3 Krugman's "Law of Motion of the Economy" 

Krugman (1992) (see also Krugman (1996b) for the results) discusses a dynamic spatial model, 
where "workers are assumed to move toward locations that offer them higher real wages" . Krugman 
makes no attempt "to model the moving decisions explicitely" , but he assumes the following "law 
of motion of the economy" : 

~jf = P X j ( u i ~ Q ) ( 9 ) 

Here, Xj is the "share of the manufactoring labor force in location j" and Uj is the real wage at 
location j. Co is the "average real wage", which is defined as: 



= 1 (io) 
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Here, the sum goes over all different regions, where j serves as a space coordinate, and it is assumed 
that "at any point in time there will be location-by-location full employment" . The assumed law of 
motion of the economy, eq. (9), then means that "workers move away from locations with below- 
average real wages and towards sites with above-average real wages" . 

In the present form, the assumed law of motion of the economy involves some shortages: 



1. A constant total number of employed workers is assumed, a change of the total number or 
unemploment is not discussed. 

2. Employed workers are considered to move. In fact, if a worker wants to move to a place which 
offers him higher wages, he first has to be a free, that means an unemployed worker, then 
moves and then has to be reemployed at the new location, again. This process is completely 
neglected (or else, it is assumed that it occurs with infinite velocity). 

3. Workers at location j always exactly know about the average wage in the system. It is not 
explained where they get the information about the average wage from. 

4. Workers move immediately if their wage is below the average, regardless of the migration 
distance to the places with higher wages, and with no doubt about their reemployment. 



In the following, Krugman's law of motion of the economy shall be derived from the general dynamic 
model presented in Sect. 2., in order to elucidate the implicit assumptions which lead to eq. (9). 
The derivation is based on four approximations: 

(i) In his paper, Krugman (1992) does not discuss unemployment. But, with respect to the dynamic 
model presented in Sect. 2, the unemployed agents exist in an overwhelming large number. Thus, the 
first assumption to derive Krugmans law of motion is, the local change of the density of unemployed 
agents due to hiring and firing can be simply neglected. 

(ii) Further, it is assumed that the spatial distribution of unemployed agents is in a quasistationary 
state. This does not mean that the distribution does not change, but that the distribution relaxes 
fast into a quasistationary equilibrium, compared to the distribution of the employed agents. 

With the assumptions (i) and (ii), eq. (8) for the density of the unemployed agents reduces to: 



d_ 

dt 



d d 2 







Integration of eq. (11) leads to the known canonical distribution: 

stat^ ^ . . - ex P [/ f(r, t)/D n dr] 



n 



(r,t) 



n 



(exp [Jf(r,t)/D n dr] 



(11) 



(12) 



exp 



J f(r, t)/D n dv^i = - jf exp f(r', t)/D n dr' 



dr 
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where the expression (...) describes the mean value, and n = N\/A is the mean density of unem- 
ployed agents. 

(iii) For a derivation of Krugmans equation, we now have to specify f(r,t) in eq. (12). Here, it is 
assumed that a single unemployed agent which migrates due to eq. (2), responds to the total local 
income in a specific way: 

f(r,t) = -^ln[u(r)l(r,t)} (13) 

Eq. (13) means that the migrant is guided by local gradients in the total income, with uj(r) being 
the local income and l(r, t) the local density of employed workers. We note here again, that the 
migrant does not count on information about the highest global income, he only "knows" about 
his vicinity. With assumption eq. (13), we can rewrite eq. (12) using the discrete notation, which 
is preferred by Krugman: 

nf*{t) = n "* l >® (11) 

EjWj ij{t) I Ej 



The corresponding equation (7) for lj(t) reads in the discrete notation: 

^l j {t) = k + n j {t) - k-l.it) (15) 

Since the Krugman equation deals with shares instead of densities, we have to divide eq. (15) by 
^Zjlj{t), which leads to: 

(iv) Krugman assumes that the total number of employed workers is constant, so we use this 
assumption to replace the relation between the hiring and the firing rate in eq. (16). Eq. (15) 
yields: 

* F = §£ < 17) 

rij in eq. (16) is now replaced by the quasistationary value, eq. (14). Inserting further eq. (17) into 
eq. (16) and using the definitions of the average wage u>, eq. (10), and n = J2j n j/J2j> we finally 
arrive at: 

±\ j (t) = ^\ j (u j -Q) (18) 

Eq. (18) is identical with Krugmans law of motion of the economy, eq. (9), if the prefactor p in eq. 
(9) is identified as: p = k~ /Co. Hence, p is a slowly varying parameter, which depends both on the 
firing rate, which determines the time scale, and on the average wage, iD, which may change in the 
course of time. 

It is an interesting question whether the assumptions (i)-(iv) which lead to Krugmans law of motion 
in the economy, have some practical evidence in an economic context. From a more theoretical 
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perspective, it is noteworthy that Krugmans equation, (9), (18) has an obvious analogy to a selection 
equation of the Eigen-Fisher type, X{ = Xi(Ei — (Ei)). Here, Ei is the fitness of species i and (Ei) is 
the mean fitness representing the global selection pressure. It can be proved (Feistel and Ebeling, 
1989) that this equation describes a competition process which finally leads to a stable state with 
only one surviving species. The competition process may occur on a very long time scale, but 
asymptotically a stable coexistence of many different species is impossible. 

In the economic context discussed by Krugman (1992), eq. (9) implies that finally all workers are 
located in one region j, where the local real wage Uj is equal to the mean real wage, u>. However, 
in his paper, Krugman (1992) (see also Krugman (1996b)) discusses computer simulations with 12 
locations (on a torus) which show the stable coexistence of two (sometimes three) centers, roughly 
evenly spaced across the torus, with exactly the same number of workers (cf. Fig. 7 in Krugman 
(1992) and Fig. 2.2. in Krugman (1996b)). The paper does not provide information about the time 
scale of the simulations, and stochastic influences (apart from a random initial configuration) are 
not considered. 

Krugmans model of spatial aggregation includes of course some more complex economic assump- 
tions, such as consideration of transportation costs, distinction between agricultural and manifac- 
tured goods, price index etc. We are not going to discuss whether the stable coexistence of two 
centers in Krugmans computer simulations might result from those specific economic assumptions 
in his model or from the lack of fluctuations (which could have revealed the instability of the (de- 
terministic) stationary state). Due to Krugman (1992), already a slight variation in the parameters 
always led to a single center. 

In the following, we will focus on the more interesting question, whether the general dynamic 
model of migration and employment introduced in Sect. 2, is able to produce a stable coexistence 
of different economic centers under the presence of fluctuations. This may allow us to overcome 
some of the shortages involved in Krugman equation. 

4 Migration and Employment of Workers due to Wage Differences 
4.1 Effective Diffusion 

The derivation of Krugmans equation was based on the assumtion (ii) that the time scale for hiring 
and firing of workers is more determining than the time scale of migration. If we explicitely consider 
unemployment in our model, this assumption implies that there are always enough unemployed 
workers which can be hired on demand. A growing economy, however, might be determined just by 
the opposite limiting case: It is important to attract workers/consumers to a certain area before 
the output of production can be increased. Hence, the time scale for the dynamics is determined by 
migration processes. That means the spatial distribution of the employed agents can be assumed 
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in a quasistationary equilibrium compared to the spatial distribution of the unemployed agents. 

J^M)=0 l stat (r,t) = ^n(r,t) (19) 

Hence, the local density of employed agents can be expressed as a function of the local density of 
the unemployed agents available, which itself changes due to migration on a slower time scale. 

We further assume that the unemployed agent who is able to migrate, eq. (2), responds to local 
gradients in the real wages (per capita), instead of gradients in the total income, as assumed for 
Krugmans equation. The migrant tries to move towards places with a higher wage, but again he 
only counts on information in the vicinity. Hence, the guiding force f(ri) is determined as follows: 

f(r,t) = ^u;(r,t) (20) 

For a further discussion, we need some assumptions about the local distribution of the real wages, 
Lo(r,t). It is reasonable to assume that the local wages may be a functional of the local density 
of employed agents, u(r, t) = uj{l(r,t)}. Some specific assumptions about this dependence will be 
discussed in the next section. With eq. (19), we can then rewrite the spatial derivative for the wages 
as follows: 

du(r, t) _ 6u>{l(r, t)} dl(r, t) _ 6u{l(r,t)} k+_ dn(r, t) 
dr 51 dr 51 k~ dr 

where 5 denotes the functional derivative. Using the equations (19), (20), (21), the Fokker-Planck 
equation for the change of the density of the unemployed agents, eq. (8) can now be rewritten as 
follows: 

£■*•*> = I Ov*^) <*> 

The r.h.s. of eq. (22) now has the form of a usual diffusion equation, with D e ff being an effective 
diffusion coefficient: 

k + 5lu 

D eff = D n - — — n(r,t) (23) 

Here, D n is the "normal" diffusion coefficient of the unemployed agents, eq. (8). The additional 
terms reflect that the unbiased diffusion is changed because of the response of the migrants to local 
differences in the wage distribution. As we see from eq. (23), there are two contradicting forces 
determining the effective diffusion coefficient: the normal diffusion, which keeps the unemployed 
agents moving, and the response to the wage gradient. 

If the local wage decreases with the number of employed agents, then 5lo/51 < 0, and the effective 
diffusion increases. That means unemployed agents migrate away from regions where employment 
may result in an effective decrease of the marginal output. However, if Slo/51 > 0, then their wage 
effectively increases in regions with a larger number of employed agents, and they are attracted 
to these regions. As we see in eq. (23), for a certain positive feedback between the local wage and 
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the employment density, the effective diffusion coefficient can be locally negative and unemployed 
agents do not leave these areas once they are there. Economically speaking, these agents stay there, 
because they may profit from the local increase in the employment density. 

We want to emphasize that this interesting dynamic behavior has been derived without any explicit 
economic assumptions. But for a further discussion of the model, the three remaining functions 
which are unspecified by now: (i) 5u/5l, (ii) k + , (hi) k~ , have to be specified, and that is of course 
where economics comes into play. 

4.2 Determination of the Production Function 

In order to determine the economic functions, we refer to a perfectly competitive industry (where 
"perfect" means complete or total). In this standard model (see e.g. Case and Fair (1992)), the 
economic system is composed of many firms, each small relative to the size of the industry. These 
firms represent seperate economies sharing common pools of labor and capital. New competitors 
can freely enter/exit the market, hence the number of production centers is not limited or fixed. 
Further, it is assumed that every firm produces the same (one) product and every firm uses only 
one variable input. Then, the maximum profit condition tells us that "firms will add inputs as long 
as the marginal revenue product of that input exceeds its market price." In the case of labor as 
variable input, the price of labor is the wage, and a profit maximizing firm will hire workers as long 
a the marginal revenue product exceeds the wage of the added worker. 

The marginal revenue product MRP = MP x P is the additional revenue a firm earns by employing 
one additional unit of input, where P is the price of output and MP is the marginal product. In 
a perfectly competitive industry, however, no single firm has any control over prices. Instead, the 
price results from the interaction of many suppliers and many demanders. In a perfectly competitive 
industry, every firm sells its output at the market equilibrium price, which is simply normalized to 
one, hence the marginal revenue product is determined by MP. 

The marginal product, MP, can be derived from a production function - Y(r,t) = A(r,t) g(r,t), 
which describes the relationship between inputs and outputs (e.i. the technology of production). 
Usually, the production function may include the effect of different inputs, such as capital, public 
goods or natural resources. Here, we concentrate only on one variable input, labor. Thus, g(r, t) is 
assumed a Cobb-Douglas production function l@(r,t) with a common (across regions) exponent (3, 
l(r,t) being the local density of employees. The exponent (3 describes how a firms output depends 
on the scale of operations. Increasing returns to scale are characterized by a lower average cost with 
an increasing scale of production, hence (3 > 1 in this case. On the other hand, if an increasing 
scale of production leads to higher average costs (or lower average output), we have the situation 
of decreasing returns to scale, with (3 < 1. In the following, we will restrict the discussion to (3 < 1, 
common to all regions. 
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The prefactor A(r, t) represents economic details of the level of productivity. We assume that these 
influences can be described by two terms: A(r, t) = A c + A u . A c = const, should summarize those 
output dependences on capital, resources etc., which are not explicitely discussed here. The new 
second term, A u , considers cooperative effects which result from interactions among the workers. 
Since all cooperative effects are non-linear effects, A u should be a non-linear functional of l(r,t): 
A u = A u {l(r,t)}. 

Hence, the production function depends explicitely only on l(r,t), now: 

Y{l(r, t)} = [A c + A u {l(r, t)}] ^(r, t) (24) 

The marginal product (MP) is the additional output by adding one more unit of input, i.e. MP = 
dY/dl, if labor is input. The wage of a potential worker will be the marginal product of labor: 

wWr , t)} = «M (25) 

If a firm faces a market wage rate of oj* (which could be bound e.g. by minimum wage laws), then, 
in accordance with the maximum profit condition, a firm will hire workers as long as: 

ol 

To complete our setup, we need to discuss how the prefactor A u depends on the density of employees, 
l(r, t). Here, we assume that the cooperative effects will have an effect only in the intermediate range 
of I. For small production centers, the synergetic effect resulting from the mutual stimulation among 
the workers is too low. On the other hand, for very large production centers, the advantages of the 
cooperative effects might be compensated by the disadvantages of the massing of agents. Thus we 
will assure, that Y ~ I 13 both in the limit I — > and I — > oo. These assumtions are concluded in the 
ansatz: 

A u {l(r,t)} ~exp{u{Z(r,t)}} (27) 
where the utility function u(l) describes the mutual stimulation among the workers in powers of I: 

u(l) = a + ail + a 2 l 2 + ... (28) 

The series will be truncated after the second order. The constants ai characterize the effect of 
cooperation, with a$ > 0. Especially, the case a\ > 0, a 2 < considers saturation effects in the 
cooperation, e.i. the advantages of cooperation will be compensated by disadvantages of crowding. 
This idea implies that there is an optimal size for taking advantages of the cooperative effect, 
which is determined by the ratio of a\ and a 2 . If one believes that cooperative effects are an always 
increasing function in I, then simply a 2 > can be assumed. 

If cooperative effects are neglected, i.e. a\ = a 2 = 0, then we should obtain the "normal" production 
function, Y(l) = Al 13 . The constant A is determined by a relation between A c and ao which results 
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from the eqs. (24), (27). Within our approach, the constant A c is not specified. If we, without 
restrictions of the general case, choose: 



A = A c + exp{a } = 2A C 



(29) 



then the production function, Y(l), eq. (24), with respect to cooperative effects can be expressed 
as follows: 



Y{l(r,t)} 



A 



1 + exp ( a\ I + Ct2 I' 



(30) 



Fig. 1 presents the production function, eq. (30) dependent on the density of employees, which can 
be compared with the "normal" production function, Y{1) = Ar . Clearly, we see an increase in 




Figure 1: Production function Y(l) vs. density of employed workers, I. The solid line results from eq. 
(30) (parameters: A = 2, a\ = 0.06, a,2 = —0.035, (3 = 0.7). The dashed line shows the production 
function without cooperative effects (ai = 02 = 0). 

the total output due to cooperation effects among the workers. If we assume a\ > 0, 02 < 0, this 
increase has a remarkable effect only in an intermediate range in I. Both for I — > and I — > 00 the 
cooperative effect vanishes. 

Once the production function is determined, we also have determined the local wage uo{l(r,t)}, eq. 
(25), as a function of the density of employees: 



u{l(r,t)} = — [l + exp(ai/ + a 2 / 2 )] pi?- 1 

A / \ 
+ - exp (01 1 + a 2 I 2 ) (ai + 2a 2 /) f (31) 

Hence, also the derivative Su /SI, used for the effective diffusion coefficient, D e ff, eq. (23) is deter- 
mined. Fig. 2 shows both functions dependent on the density of employees. 
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Figure 2: (left) Wage to, eq. (31), (right) derivative 5u)/5l vs. density of employed workers, I. The 
solid lines show the functions with, the dashed lines without cooperative effects (for the parameters 
see Fig. 1). 

The left part of Fig. 2 indicates that, within a certain range of /, the derivative Su/Sl could be 
indeed positive, due to cooperative effects. With respect to the discussion in Sect. 4.1., this means 
that the effective diffusion coefficient, eq. (23), can be possibly negative, i.e. unemployed workers 
will stay in these regions. 

4.3 Determination of the Transition Rates 

The transition rate k + for hiring unemployed workers is implecitely already given by the conclusion, 
that firms hire workers as long as the marginal revenue product exceeds the wage of the worker. In 
accordance with eq. (26) we define: 

fc+ = k + {l(r, t)} = a exp j — t)} - a;*} (32) 

Here, a determines the time scale of the transitions. k + , which is now a function of the local 
economic situation, is significant larger than the level of random transitions, represented by a, only 
if the maximum profit condition allows hiring. Otherwise, the hiring rate tends to zero. 

The firing rate k" could be simply determined opposite to k + . Then, from the perspective of the 
employee, firing is caused by the local situation of the economy, i.e. by external reasons: workers 
are fired if SY/51 < uj* . A more refined description, however, should consider also internal reasons: 
a worker cannot only loose his job, he can also quit his job himself for better opportunities, e.g. 
because he wants to move to a place where he earns a higher wage. 
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It is reasonable to assume that the internal reasons depend again on spatial gradients in the wage. 
Due to eq. (2) the unemployed agent migrates while guided by gradients in the wage. The employed 
agent at the same location may have the same information. Noteworthy again, this is only a local 
information about differences in the wage. If the local gradients in the wage are small, the internal 
reasons to quit the job vanish, and the firing depends entirely on the (external) economic situation. 
However, if these differences are large, the employee might quit his job for better chances. 

We note that the latter process was already considered in Krugmans law of motion in the economy, 
Sect. 3. Different from the assumptions involved in Krugmans eq. (9), here the process: employment 
— ► unemployment — > migration — ► reemployment is explicitely modeled. It does not occur with 
infinite velocity, and there is no guarantee for reemployment. 

Hence, we define the "firing rate" which describes the transition from an employed to an unemployed 
agent, as follows: 



k =k {l(r,t)} =a expj- ^ ^ - c^*j + q-^rj 



(33) 



The additional parameter q can be used to weight the influence of spatial gradients on the employee. 

Eventually, in this section we have determined the variables f(r,t), k + and k~ via a production 
function Y{l(r,t)}, which represents certain economic assumptions. Now, we can turn back to the 
dynamic model described in Sect. 2, which now can be solved by means of computer simulations. 



5 Numerical Simulations 

5.1 Stochastic Simulation Technique 

Before presenting the results of the computer simulations, the simulation technique should be 
shortly discussed. The computer program has to deal with three different processes, which have to 
be discretized in time for the simulation: (i) the movement of active Brownian particles with 9i = 1 
due to the overdamped Langevin eq. (2), (ii) the transition of the particles due to the rates, k + , 
eq. (32), k~ , eq. (33), and (iii) the generation of the field u(r,t). 

Considering eq. (2), the new x-position of a particle % with B\ = 1 on the two dimensional surface 
at time t + At is given by: 

At + ^2D n At GRND (34) 

The equation for the y-position reads accordingly. At is the non-constant time step, which is 
calculated below. D n = e is the diffusion coefficient, and GRND is a Gaussian random number 
with mean equals zero and standard deviation equals unity. 



Xi (t + At) =Xi(t) + — 
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In oder to calculate the spatial gradient of the wage field u, eq. (31), we have to consider its 
dependence on l(r, t), which is a local density. The density of employed agents is calculated assuming 
that the surface is divided into boxes with the spatial (discrete) indices u, v and unit length As. 
Then, the local density is given by the number of agents with Oi = inside a box of size (As) 2 : 

1 rx u +As ry v +As 

l(r,t) ^ l uv (t) = —— / dx'dy' C (x',y') (35) 

(AS) 2 J Xu Jy v 

We note that As » Ax, with Ax = x(t + At) — x(t) being the spatial move of the migrating agent 
into x-direction during time step At. That means that the migration process is really simulated as 
a motion of the agents on a two dimensional plane, rather than a hopping process between boxes. 

Using the box coordinates u, v, the production function and the wage field can be rewritten as 
follows: 

Y {Z(r,t)} Y{l uv {t)} = Y uv (t) ; w(r) u uv = (36) 
The spatial gradient is then defined as: 

diOuv _ Wu—l,v ~ &u+l,v _ dto uv _ UJ UjV -i — u UtV+ i 

~dx~ ~ 2A~s ' ~dy~ ~ 2A~s 1 ' 

where the indices u—1, u + 1, v — 1, v + 1 refer to the left, right, lower and upper boxes adjacent to 
box uv. We further note that for the simulations periodic boundary conditions are used, therefore 
the neighboring box is always specified. 

Using the discretized versions, eq. (35), (36), (37), the transition rates can be reformulated as k£ v , 
k~ v accordingly. They determine the average number of transitions of a particle in the internal state 
located in box uv, during the next time step, At. In a stochastic simulation however, the actual 
number of reactions is a stochastic variable, and we have to assure that the stochastic number of 
reactions (i) does not exceed the actual number of particles available during At, and (ii) is equal 
to the average number of reactions in the limit t — > oo. 

This problem can be solved by using the stochastic simulation technique for reactions, which defines 
the appropriate time step, At, as a random variable. Let us assume that we have exactly Aq particles 
with the internal state 0j = 1 in the system at time t = to- Then the probability P(Aq,t ) equals 
one, and the probability for any other number N[ is zero. With this intial condition, the master 
equation to change Aq reads: 

^Qt^ =J2(kt + K)P(N 1 ,t) ; P(N 1 ,to) = l (38) 
i=i 

where ./V is the total number of particles, and kf , k~ are the transition rates for each particle. The 
solution of this equation yields: 



P(N u t-to) ~ exp(-LJ^) 
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1 N Ni N 

I™ = — J E^ = E^ + + E K (39) 
J2 hi i=1 i=1 i=JVi+l 

1=1 

Here, i m is the mean life time of the state N\. For i — to <C i m , the probability P(Ni,t) to find still 
iVi is almost one, but for t — to ^ t m this probability goes to zero. The time when the change of 
Ni occurs, is most likely about the mean life time, t — to ~ t m . In a stochastic process, however, 
this time varies, hence, the real life time t — to = t is a randomly distributed variable. Since we 
know that P(N\,t) has values between [0, 1], we find from eq. (39): 

r = t-to = -tm ln{i?iVD[0,l]} (40) 

(Zero actually has to be excluded.). That means, after the real life time r, one of the possible 
processes which change N\ occurs. Each of these processes has the probability: 

N 

i=i 

N 

p(9 i = l,t^9 i = 0,t + r) = k+/J2h (41) 

i=i 

Thus, with a second random number RND[0,J2 hi] it will be determined which of these possible 
processes occurs. It will be the process q which satisfies the condition: 

q-l q 

J2h< RND ; k i > RND ( 42 ) 
i=i i=i 

For the transition probabilities to change iVi, we find in particular: 

N 

Pr6b[N u t^ iVi + M + r] = E ^ rt ™ 

i=i 

Ni 

Prob[N u t ^m-^t + T] = E^ +t ™ ( 43 ) 

1=1 

Obviously, the sum over these probabilities is one, which means, during the time intervall r one of 
these processes occurs with certainity. Further, using the definition of r, we see that 

* = E4 = I>' ( 44 ) 
i=i i=i 

yields asymptotically, i.e. after a large number of simulation steps, z S> 1. 

Hence, determining the time step as At = r, eq. (40), ensures that only one transition occurs during 
one time step and the number of transitions does not get out of control. Further, in the asympotitc 
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limit, the actual (stochastic) number of transitions is equal the average number. The (numerical) 
disadvantage might be that the time step At is not constant, so it has to be recalculated after each 
cycle before moving the particles with respect to eq. (34). This may slow down the speed of the 
simulations considerably. 

Let us conclude the procedure to simulate the movement and the transition of the particles: 

1. calculate the density of particles l uv (t), eq. (35) 

2. calculate the production function Y uv {t) and the wage field w uv (t), eq. (36) 

3. calculate the sum over all possible transitions, which determines the mean life time t m , eq. 
(39), of the current system state 

4. calculate from the mean life time the actual time step At = r, eq. (40) by drawing a random 
number 

5. move all particles with Bi = 1 according to the Langevin equation (34) using the time step r 

6. calculate which one of the particles undergoes a transition by drawing a second random 
number, eq. (42) 

7. update the system time: t = t + r and continue with 1. 

5.2 Computer Simulations of Spatial Economic Aggregation 

With the determination of the three functions f(r,t), k + and k~ with respect to some economic 
assumptions, we have completed our dynamic model, described in Sects. 2 and 4. In this section, we 
want to discuss some features of the dynamics by means of a computer with 500 active Brownian 
particles. Initially, every particle is randomly assigned an internal parameter (either or 1), and 
a position on the surface, which has been divided into 10 x 10 boxes of unit length As = 1. The 
diffusion coefficient, which describes the mobility of the migrants, is set to 0.01 (in arbitrary units), 
So, a simple Brownian particle would approximately need a time of t = 50 for a mean spatial 
displacement of 1 (which is the spatial extension of a box). The minimum wage to* is set to 0.1. 
Further, a = 1 and 9 = 1, for the remaining parameters see Fig. 1. 

In the following we will restrict the discussion to the spatio-temporal evolution of the densities of 
employed and unemployed agents. Other quantities of interest, such as the spatio-temporal wage 
distribution, the production function and the local values of the "hiring" and "firing" rates have 
been of course calculated for the simulation, but will be discussed in a subsequent paper, which 
also investigates different parameter sets. 

Fig. 3 shows snapshots of the evolution of the spatial density of employed agents, while Fig. 4 shows 
the corresponding spatial density of unemployed agents. Fig. 5 presents the evolution of the total 
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number of employed and unemployed agents, in terms of the total share xg = N$/N with respect 
to eq. (6). 





Figure 3: Spatial density of employed agents, l(u,v), for different times t (for the parameters and 
initial conditions of the simulation, see text). 

The simulation shows that the evolution of the densities occurs on two stages: 

(i) During the first stage (for t < 1.000) we have a significant higher share of employed agents (up 
to 70 percent of the agents community) . They are broadly distributed in numerous small economic 
centers, which have an employment density of about 5. This can be understood using Figs. 1, 2. 
The increase of the output due to cooperative effects allows the fast establishment of numerous 
small firms, which count on the mutual stimulations among their workers and in the beginning offer 
higher wages (compared to the model without cooperative effects). 

This growth strategy, however, is not sufficient for output on larger scales, because the marginal 
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Figure 4: Spatial density of unemployed agents, n(u,v), for different times t (for the parameters 
and initial conditions of the simulation, see text). 



product drastically decreases (before it may increase again for an above critical employment den- 
sity). As a result the wages may fall (possibly below the minimum wage rate, u*), which prevents 
further growth. So, the first stage is characterized by the coexistence of numerous small economic 
centers. 

(ii) During the second and much longer stage (for t > 1.000), some of these small centers overcome 
this economic bottleneck, which allows them to grow further. In the simulations, this crossover 
occured between t = 500 and t = 800. In the model considered here, we have only one variable 
input, labor, so the crossover is mainly due to stochastic fluctuations in Z(r, t). 

After the bottleneck, the marginal output increases again, and also the wage increases with the 
density of employment (cf. Figs. 1, 2). This in turn effects the migration of the unemployed agents, 
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Figure 5: Total share x$ of employed agents (□) and unemployed agents (O) vs. time t. The marks 
indicate simulation data which correspond to Figs. 3, 4. 

which, due to eq. (2), is determined by two forces: the local attraction to areas with higher wages, 
and the stochastic influences which keep them moving. With du/dl > 0, the attraction may exceed 
the stochastic forces, so the unemployed agents are bound to these regions, once they got there. 
This, in turn is important for the further growth of these economic centers, which needs agents to 
hire. As a result, we observe on the spatial level the concentration of employed and unemployed 
agents in the same regions. 

For their further growth, the economic centers which overcome the bottleneck, locally attract the 
labor force at the expense of the former small economic centers. As a result of the competition 
process, these small centers, which previously employed about 70 percent of the agents, dissapear. 
In the new (larger) centers however, only 60 percent of the agents can be employed (cf. Fig. 5), 
so the competition and the resulting large-scale production effectively results also in an increase 
of unemployment. But the employment rate of 60 percent seems to be a stable value for the given 
parameter setup, since it is kept also in the long run, with certain fluctuations. 

The concentration process described above, leads to the existence of different extended major eco- 
nomic regions, shown in Figs. 3, 4, which each consist of some subregions (in terms of boxes). In 
the long run (up to t=100.000) we find some remarkable features of the dynamics: 

(a) a stable coexistence of the major economic regions, which keep a certain distance from each 
other. This critical distance - which is a self- organized phenomenon - prevents the regions from 
detracting each other. In fact, finally there is no force between these regions which would pull 
off the employed or unemployed agents and guide them over long distances to other regions. 
Thus, the coexistence is really a stable one. 
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(b) a quasi- stationary non- equilibrium within the major economic centers. As we see, even in 
the long run, the local densities of employed and unemployed agents do not reach a fixed 
stationary value. Within the major regions, there are still exchange processes between the 
participating boxes, hiring and firing, attraction and migration. Hence, the total share of 
employed agents continues to fluctuate around the mean value of 60 percent (cf. Fig. 5). 

6 Conclusions 

In this paper, we have proposed a simple dynamic model to descibe migration and economic aggre- 
gation within the framework of Active Brownian particles. We consider two types of particles (or 
agents); employed agents, which are immobile, but (as the result of their work) generate a wage 
field, and the unemployed agents which migrate by responding to local gradients in the wage field. 
Further a transition between employed and unemployed agents (hiring and firing) is considered. 

The economic assumptions to be used in the model, are concluded in only three functions: (a) the 
derivative doj/dl, which describes how the local wage depends on the employment density, (b) k + , 
(c) k~ , which are the hiring and firing rates. These functions can be determined using a production 
function, Y(l), which describes the output of production dependent on the variable input of labor. 
Our ansatz for Y(l) specifically counts on the influence of cooperative interactions between the 
agents on a certain scale of production. 

As the result of our dynamic model, we find the establishment of distinct economic centers out 
of a random initial distribution. The evolution of these centers occurs in two different stages. 
During the first stage, small economic centers are formed based on the positive feedback of mutual 
stimulation/cooperation among the agents. During the second and much longer stage, some of 
these small centers grow at the expense of others, which leads to the interesting situation that in 
our model economic growth and decline occur at the same time, but at different locations, which 
results in specific spatial-temporal patterns. The competition finally leads to the concentration of 
the labor force in different extended economic regions. This crossover to large-scale production is 
accompanied by an increase in the unemployment rate. 

Although the extended economic regions are in an internal, quasistationary non-equilibrium, we 
observe the stable coexistence between these regions. This is an important result, which agrees 
e.g. with the central place theory of economics (Christaller, 1933, Losch, 1940). Different from an 
attempt by Krugman (1992), who also focused on this problem, we find that (i) the coexistence 
is stable, even at the presence of fluctuations, (ii) the centers not necessarily have to have the 
same number of employed agents, to coexist, (iii) the dynamics does not simply converges into an 
equilibrium state, but the centers exist in an quasistationary non-equilibrium state and still follow 
a stochastic eigendynamics. 
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One question, which will be tackled in a forthcoming paper, is about the influence of the different 
parameters. We have seen in the simulations, that cooperative effects are able to initiate an increase 
in the economic output, which leads to the establishment of many small firms. But the final situation 
does not include any of those small firms and seems to be independent of this intermediate state. 
This is mainly due to the fact, that we considered only one product. A more complex production 
function with different outputs, however, should also result in the coexistence of small (innovative) 
firms and large scale production. 

Another important issue is to understand the relation between the response of the agents to local 
gradients in the wage field, on one hand, and stochastic influences, on the other hand. Most likely, 
there exist critical parameters which may describe a "phase transition" within the agents behavior, 
which eventually results in different types of spatial-temporal coexistence. 

Finally, one may argue that the simple model proposed here, is not the simplest one to describe 
the economic aggregation. Hence, the minimalistic agents discussed in Sect. 1, could be even more 
minimalistic. One example is the response of the employed agents to gradients in the wage field, 
expressed in the parameter q which enters the firing rate, k~. q = simply means that internal 
reasons are neglected, and the agent is fired only due to external reasons. This might be also a 
reasonable assumption, however, it includes a pitfall. In a growing economy, eventually almost all 
unemployed agents might be hired at some locations. But, as long as dY/dl > gj*, employed agents 
are almost never fired, thus there is a shortage of free agents for further growth. Eventually the 
dynamics sticks in a dead-lock situation, because agents cannot move even if there are locations 
in their neighborhood which may offer them a higher wage. Therefore, it is reasonable to choose 
q > 0, but it is important to understand the role of q, for instance, for the final patterns observed. 

So, we conclude that our model may serve as a toy model for simulating the influence of differ- 
ent social and economic assumptions on the spatial-temporal patterns in migration and economic 
aggregation. 
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